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TWO INTERFERING QUEUES IN PACKET-RADIO NETWORKS 
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ABSTRACT 


We consider several classes of interfering queues that appear in packet-radio 
networks. We analyse the class of systems where one of the queues is given full 
priority and obtain an expression for the joint probability distribution of the queue 
lengths. For ALOHA-tvpe systems with two symmetric queues we calculate the average 
packet waiting time and queue lengths and for symmetric systems with an arbitrary 
number of subscribers we develop a method to approximate these quantities. The 
approximation turns out to be close to the analysis and simulation results. 
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1. Introduction 


The present study was motivated by the problem of investigating the behavior 
of random-multiple-access systems and of packet-radio networks. These systems are 
characterized by the fact that a number of radio stations exchange digital informa¬ 
tion by using a distributed random access algorithm on a common radio channel. In 
such situations, whenever a given station attempts transmission of a packet to 
another station, the attempt may be unsuccessful, in which case the packet must be 
retransmitted. In addition to channel noise, unsuccessful transmissions occur 
because of interference from another station trying to send a packet over the 
common channel at the same time or by the fact that the intended receiver is itself 
in transmission mode, in which case it is not able to detect incoming packets. The 
fact that the activity at one node affects the behavior of the queue at other nearby 
stations gives rise to statistical dependence between the queues at the nodes. 

Unfortunately, in general,the queue length statistical dependence is quite 
complicated and there is little hope to obtain explicit analytical results for general 
topology networks. The purpose of this paper is to present several analytic as well 
as approximation results for certain classes of interfering queues. We assume 
throughout the paper that all packets have equal length and that the time is divided 
in slots corresponding to the transmission time of a packet. A station may start 
packet transmissions only at the beginning of a slot and the distances between stations 
are assumed to be such that propagation delay is negligible. Also, we neglect channel 
noise and assume no channel errors. 

Because of the difficulty in the analysis of dependent queues of the type 
introduced above, even the case of two queues cannot be treated analytically in the 
general case. However, we consider here two classes of systems with two dependent 
queues where such results can be obtained. The first is the case when the length of 
one of the two queues is not allowed to decrease, unless the other queue is empty. A 
two node ALOHA system where one of the nodes is given full priority is an example for 






such a situation; other examples are given in Sec. 2. No other restrictions are 
necessary in order to allow for analytical solution of this class; in particular, 
the inputs may have arbitrary distributions and need not be independent processes. 

For this class of problems we present a general methoa for deriving the generating 
functions of the queue lengths and of the average delay times. Then these general 
results are applied to three special cases of packet-radio networks that can be shown 
to belong to the considered class of systems. 

The second class of problems for which we can obtain explicit analytical 
results is the case of a two node symmetric ALOHA system. For this situation we 
cannot obtain the queue length probability distribution (or generating function), 
but we give a method for calculating the average queue length and hence the average 
time delay. The results are given in Sec. 3. 

Since, as said before, explicit analytical results are hard to obtain for more 
general situations, another way to approach the problem is to obtain good approxima¬ 
tions. An approximation method applicable to a symmetric ALOHA-type system with 
arbitrary number of stations is introduced in Sec. 4. We obtain there the approximat 
average queueing delay in such systems and compare this with the exact result for two 
node networks and with simulation results for larger networks. 

Discrete time systems involving interfering queues have been rarely treated in 
the literature. In [2] a loop system using Asynchronous Time-Division Multiplexing 
has been analysed. In this system user i may use a slot for transmission only if 
all users l,2,...,i-l have nothing to transmit. In [3], the author investigates 
the case of two queues in tandem where each queue always attempts transmission 
provided it has a message in the buffer. Both systems considered in [2] (for two 
users) and [3] are cases that belong to the first general class of systems considered 
in the present paper. In [4] another system of two interfering queues is considered, 
whereby only one event (i.e. an arrival or departure) may occur during a given slot. 
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This system has been shown to have the product form solution. Finally, we may 
mention [5] and [6], where a slotted ALOHA network with finite number of users has 
been examined and a method was suggested for obtaining an approximate solution for 



mis svstem. 
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2. Analytical Results 

In this section we consider a class of discrete-time queueing systems consisting 
of two queues with the following properties : packets arrive randomly at the queues 
from two sources, that in general may be correlated. Let A^(t) and A 2 Ct) be the 
number of packets entering node 1 and node 2 from their corresponding sources in 
the time interval (t, t+1]. The input process [(A (t), A 2 Ct)] is assumed to be a 
sequence of independent identically distributed random vectors with integer-valued 
elements. Let 

CO 00 

a(i,j) = Prob(A (t)=i,A ? (t)=j) ; £ £ a(i,j) = 1 . (1) 

1 i=0 j=0 

and 

F(x,y) = E[x Al(t) y A2(t) ] » I [ a(i,j)x i y i (2) 

i=0 j«0 

We assume that F(x,v) cannot be x-independent, namely that packets arrive at the 
first queue with nonzero probability and that the queues have infinite buffers. 

Next, we describe the departure processes. It is assumed that no more than 
one packet may leave each queue in any given time slot and the combined departure 
process is taken to be as follows : When both queues are empty, no departures may 
occur (packets arriving during a given slot may depart only in the next one). When 
only one of the queues is nonempty, a departure from that queue may occur and the 
packet may be transferred either to the outside of the system or to the other queue. 

We denote by p^, p^ the respective probabilities when the nonempty queue is queue 
2 and by P^q. pJq the corresponding probabilities when the nonempty queue is queue 1. 
Clearly 1 - p^ - p^ is the probability that no departure occurs from the nonempty 
queue 2. The specific class of dependent queues considered in this section is 
characterized by an assumption on departures when both queues are nonempty. For this 
case it is assumed that a departure may occur only from queue 2. We denote by p^, 
pjj the probabilities that the departing packet leaves the system or joins the other 
queue (queue 1) respectively. 


i 







Consider the steady-state joint generating function of the queue lengths : 


r L.ct) L 2 (tf 

C(x,y) = lim Ex y (3) 

t-*-oo I— - ■ 

where L^Ct), L 2 (t) are the queue lengths at time t at nodes 1 and 2 respectively, 
and where we assume that the Markov chain [L^(t), L^Ct)] is ergodic, namely 
G( i 0,0)>0. For the class under consideration we can compute the function G(x,y) and 
in the Appendix it is shown that G(x,y) has the following form : 


rev v) - FTy vi • btXyjGCx.O) + c(x,y)G(0,y) + d(x,y)G(0,0) 
,>J ~ 1 ,>j x• e(x,y] 


where the functions b(x,y), c(x,y), d(x,y), e(x,y), G(y,0), G(0,y) and the constant 
G(0,0) are defined in the Appendix. 

This general form can be ma<Te more explicit for certain special cases. We 
next consider several examples of two node packet networks, where it turns out that 
the general assumptions, given earlier, characterizing the class of two dependent 
queues indeed hold. The networks under consideration are given in Fig. 1. In all 
cases the nodes share a common radio channel and are equipped with radio transmitter 
devices and in systems 2 and 3, node 2 has also a receiving device. Node 2 can either 
transmit or receive, but not simultaneously. The circle in Fig. 1 represents a station 
equipped with a radio receiver that receives packets correctly provided that there is 
no interference. Finally, instantaneous feedback to the transmitter is assumed i 
meaning that a transmitter knows at the end of the slot if the packet has been received 
correctly. 

In all three systems of Fig. 1 node 2 is assumed to have full access capability 
to the common channel. This means that it always transmit a packet when its buffer is 
not empty, while if its buffer is empty the node does not transmit and in systems 2 
and 3 it is able to receive packets transmitted by node 1. Node 1 has only partial 
access capability to the channel and its transmission policy is randomized as follows : 
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At the beginning of each slot for which its own buffer is nonempty, node 1 
tosses a coin with probability of success p, independently of any other event 
in the system and in case of success the node attempts to transmit the packet at 
the head of the queue. Both nodes are able to detect at the end of the slot if 
their transmissions were successful. At any node, if the transmission is not 
successful, either because the packet was sent to the other node while that node was 
not ready to receive it, or because of interference with a packet transmitted by the 
other node, the transmitter repeats the procedure described above. 

Since node 2 has full access capability to the channel and node 2 cannot 
receive and transmit packets at the same time, it is clear that in all cases no 
packets can leave node 1 whenever the queue at node 2 is nonempty and therefore all 


cases of Fig. 1 belong to the class of queues considered earlier in this section. 

k 

We now turn to calculate the parameters (p„, Osi.j.ksl, i+j>0} in each of 
these three systems. System 1 depicted in Fig. 1(a) represents a two node non- 


symmetric ALOHA network, where both nodes send their packets to the station. Since 
no packets are sent from one node to the other we have = pj^ = 0. When 

one of the nodes has packets to transmit while the other is empty, any attempted 


transmission is successful. Since node 2 transmits with probability 1 and node 1 
with probability p, we have p^ Q = p; p^ = 1. When both nodes have nonempty 
queues, successful transmission occurs at node 2 whenever node 1 does not attempt 
transmission and therefore 1 p^ = p. System 2, depicted in Fig. 1(b) represents a 
situation of two tandem nodes where the station that is the "sink" for the packets 


transmitted by node 2, is out of the transmission range of node 1. Therefore node 1 


cannot interfere with the transmissions of node 2. However node 2 does interfere 


with the transmissions of node 1 since when it is transmitting, it does not accept 
packets transmitted by node 1. Consequently pj^ = p; p^ = p^ = 1; 

= ^01 = Pli = System 3, depicted in Fig. 1(c) differs from system 2 only in 

that the station is in the transmission range of node 1, therefore node 1 does 
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interfere with the transmissions of node 2 in this case, and therefore the parameters 


ar: 


i 0.0-0 1 

P 10 “ p ’ P 01 = 1> P 11 " p ’ p I0 = p 01 


'll 


= 0. 


Numerical Results 

Although the results of this section hold for general input processes, the 
equations become much simpler when one considers independent Bernoulli processes. 
In this case we have 


F(x,y) = (xr x + r^ (yr 2 + r 2 ) (5) 

where r^, r 0 are the input rates. For this example we calculate for each of the 
three systems the average delays (in units of slots) T , T 2 at nodes 1 and 2 
respectively and the total average delay T in the network. This is done by first 
calculating the average queue lengths at the nodes and then applying Little's Theorem 
[7], After straightforward but tedious algebra the following results are obtained : 
System 1 (Fig. 1(a)) 

(P) 2+ r 2 P r i r 2 PP 

Tj = 1 + --—r + - j- 1 - C6) 

p(p-r 2 )-r 1 p (p-r 2 ) [p(p-r 2 )-r^] 


T 2 = ! + 


‘l p 


and 


(P-r 2 )“ 
r l T l + r 2 T 2 


r l +r 2 


(7) 


( 8 ) 


where these equations hold for p(p-r 2 ) > r^p which is the ergodicity condition in 
this system. 

In Fig. 2,3,4 we plot T^,T 2 and T respectively versus p, the transmission 
probability at node 1, for r^ =0.1 and r 2 ranges from 0.01 to 0.4. As expected 
the average delay at node 2 increases when p increases since then its transmissions 
are more frequently interfered with transmissions from node 1. More interesting is 
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the behavior of the average delay at node 1. Here p has some value for which 

is minimized (for given r^ and r 7 ) 4 When p either increases or decreases from this 

value increases. The reason is that when p becomes small, node 1 attempts 

to transmit relatively rarely, so its queue increases. When p becomes large, then 

node 1 attempts to transmit more frequently, thus interfering with the transmissions 

of node 2, and the queue lengths at both nodes are large. As we see from Fig. 4 the 

parameter p is a very critical design parameter of this system and for given values 

for r^ and r,, there exists an optimal p, denoted by p*, that minimizes the 

total average delay in the network. In Fig. 5 p* is plotted versus r^ for various 

values of r 7 . Notice that p* is much less sensitive to changes in r^ than to 

changes in r 7 . When r 7 is small, then p* should be large for all values of r , 

since interference between transmissions is rare. When r 2 increases, p* should 

decrease in order to reduce the interference. Finally, in Fig. 6, T . the minimum 

mm - 

total average delay is plotted versus y the total throughput of the system, when 
r^ = r 7 = r (clearly y = 2r). 

System 2 (Fig. 1(b)) 

The average delays are : 


T 1 = 1 + 


r l P+f 2 (1 " p? 2^ 

r 2 [p(l-r 1 -r 2 )-r i ; 


T 7 = —-— 
2 r 1+ r 2 


r T + I- 

2 1-r¬ 


and 


T = —i— T. + T, 
r l +r 2 1 2 


(9) 


( 10 ) 


(ID 


where these equations hold for p(l-r 1 -r 2 ) > r^ which is the ergodicity condition 
for this system. 

In this system node 1 does not interfere with the transmissions of node 2. 
Therefore, it is optimal to always attempt transmission at node 1 as well, namely to , 

,Cf 

take p - 1. Eq. (9) indeed shows that T^ is monotonically decreasing when p y 
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increases, and it achieves its minimum value for p ■ 1. In Fig. 7,8 and T are 

plotted respectively versus r^ for p = 1 when r 2 ranges from 0.01 to 0 . 8 . T 2 
is not plotted here since it takes values from 1 to 1.5 only when 0.01^r 2 s0.8, 

0 < r, < 0.5. From Fig. 7 we see, as expected, that Tj increases when the arrival 
rate r^ to node 1 increases. When r 0 increases, then the time spent by node 2 
in the transmitting mode increases. Therefore, packets transmitted by node 1 are 
rarely received by node 2 and the average time delay T at node 1 , increases in 
this case as seen in Fig. 7. 

System 5 (Fig. 1(c)) 

pir.+r+ (p-r ? )(p ) 2 pt^+r^r^r^) ] 

T. = 1 * --——---- 5 -I 12 ) 

[p(p-r,)- ri ](p-r 2 ) [? 1 f 2 )“-p(l-r 1 r 2 )](p ; r 2 ) 

f 

1 | V r 2 f 2 pCr 1 ^r 2 (r 1 + r 1 r 2 )][p(p-r 2 )-r 1 ] ^ 

2 r l +r - p-r 2 (p-r 2 ) [r 1 (r 2 ) 2 -p(l - r^)] 


and 
T = 


r, 

i. 

-T, 

r l +r 2 1 


To 


(14) 


where these equations hold for p(p-r 2 ) - r^ > 0 which is the ergodicity condition 
for this system. Similarily to System 1, T^, T 2 and T are plotted in Fig. 9,10,11 
for System 5. The behavior of the average delay in this system is very similar to 
that of System 1, therefore we will not give here more detailed explanations. 
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5. Symmetric Two-Node ALOHA Network Delay Analysis 

In this section a two node symmetric ALOHA network is considered. This network 
is similar to System 1 with the following modifications. Here node 2 uses the same 
channel access scheme as node 1, i.e. at the beginning of each slot, if its buffer is 
not emptVj node 2 tosses a coin, independently from any other event in the system, 
with probability of success p. According to the outcome, the node either transmits 
or remains silent during the current slot. In addition, it is assumed that the 
arrival processes to the nodes are independent Bernoulli processes with equal rates, 
denoted by0<r<l/4. Therefore we have for this system F(x,y) = (xr+r)(vr+r). By 
the same method used in the Appendix it is easy to see that for the current system 
we have : 


G(x,y)=F(x y ) p [ x l-'y- 1 )-p(2xy-x-y)lG(x,. Q ) + [yCx-l)-p(2xy-x-y)]G( 0 ,y)^pf2xy-x-y]G(0,0) 

xv-F(x,y)[(x+y)pp+xy(p'+(p)")] 

In this case we cannot obtain an explicit form for G(x,0), G(0,y), G(0,0) and hence 
for G(x,y). However we can exploit the symmetry to obtain an expression for the 
average delay in the system. If we denote by Gj(x,y), G 9 (x,y) the derivative of 
G(x,y) with respect to x and y respectively, we clearly have G^(l,l) = G-,(l,l) 
and G^(1,0) = G o (0,l). Then from (15) we obtain 


_ -> 


Gjd.l) = r 


r[p + (p)~] - P G^(1,0) 


and 


PP - r 


^[G( x ,x)] 


2 r + 


Gj (1,0)p(l-2p) r“ + 2r-4rpp 


x=l 


pp - r 


2 (pp - r) 


(16) 


(17) 


for pp > r. 

Now, if we use the fact that 


^[G(x,x)] 


x=l 


= G 1 (l,l) + G : (l,l) = 2G 1 (1,1) 


(18) 


(15) 






we can solve for G^(1,0) and hence for G^(l,l) and obtain 

,-,2 1 

(P) + 7 rp 

GjCl.l) = G 2 (l,l) = r + r----- (19) 

PP * r 

Therefore, applying Little's Theorem we obtain the average system delay 
G, (1,1) (p ) 2 + j rp 

T = - = 1 + - for pp > r . (20) 

pp - r 

2 1/2 

From (20) it is found that p* = l-(0.5r + [0.5r(l-r+0.Sr )] }/(l-0.5r) minimizes 

T for 0<r<0.25. In Fig. 6 the minimum total average delay is plotted versus 

Y, the total throughput of this system. Comparing the curves in Fig. 6 it is clear 
that the non-symmetric access scheme used in system 1 provides very slight better 
performance than the symmetric access scheme, when the arrival rates into the nodes 
are equal (the difference in the minimum delay is less than 3.5% in the range 0<y < 0.5). 
However the former scheme is unfair, giving priority to queue 2, although the arrival 
rates into the two queues are the same. 

Finally we may mention that the method presented in this section for calculating 
average delay times without obtaining explicitly the generating functions, can be used 
in many other symmetric (i.e. G(x,y) = G(y,x)) two node systems. Specifically, we can 
easily obtain average delay time in the symmetric two nodes ALOHA network for general 
arrival processes into the nodes for which F(x,y) = F(y,x). 
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4. Symmetric ALOHA Network with M Nodes : Approximate delay 

Consider an ALOHA network consisting of M nodes, that share a common radio 
channel. Assume that all nodes in the network use the same channel access scheme 
as the two nodes in the previous section, i.e. each node tosses independently a 
coin with the same probability of success p at the beginning of each slot when 
its queue is nonempty. Also assume that the arrival process at each node is a 
Bernoulli process, independent from node to node, with rate r. For this case 
we cannot obtain the exact average delay for M>2, and therefore we must consider 
approximations. The approximation method proposed here is the following. Each node 
i in the network may at any time be in one of M possible situations. Situation 
j refers to the case where j-1 nodes other than i have nonempty queues, while 
the other M-j nodes (not including i) are empty. The approximation considered here 
consists of assuming that in steady state, while in situation j, node i behaves 
as a discrete M/M/1 queue with arrival parameter r and departure parameter 
p(l-p) J_1 . If the transitions between the various situations are neglected then the 
average number of packets at node i denoted by L^ is : 


L. = >9. - v . — 

1 j = l J pC1-p)^' 1 -r 

for i=l,2,...,M and r<p(l-p) M ~\ where 
situation j. We approximate 9^ as 


0 . 
J 


( 21 ) 


is the probability of being in 


9 . 
J 



M-j r j-1 
P 


j = 1,2,. ...M 


( 22 ) 


where — approximates the probability that a node has packets ready for transmission 
and we assume independence between the nodes. Using (21) and (22) and applying 


Little's Theorem we find that the approximate total average delay T is given bv : 

a P 


T 

ap 


M 

I 

j = l 


P 


1 (1 - I)»H 

P 


1-r 

p(l-p)^ _1 -r 


(23) 
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The general formula (23) can be specialized to the case M=2, and for this case we 
can compare the approximate and exact results. For M=2 we have from (20) 


and from (23) 


T 


analysis 


,-.2 1 
(P) + 2 ^ 

PP - r 


T 

ap 


+ 


(P) 2 + rp 
PP - r 


(24) 


(25) 


T and T , are plotted versus p in Fig. 12, for various values of r. 

ap analysis r 

It can be seen that although we have used a simple approximation, it is quite close 
to the exact values. We also compare this approximation versus simulation results, 
for networks having three and four nodes. This comparison is plotted in Fig, 13 
and 14 and shows again good results. 
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Footnote; 

1. We use the notation v = 1-v. 


The subscripts 1 and 2 correspond to the derivative of the function with 
respect to the first and second variables respectively. 






In order to describe the behavior of the queue lengths, we need several 
definitions. Let D^^t) be a binary-valued random variable that takes value 
I if L (t) = 0, L 2 (t) > 0 (L^(t) and L 2 (t) are the queue lengths at time 

t at nodes 1 and 2 respectively), a departure occurs (from queue 2) and 
the packet leaves the system. Similarly D^(t) takes value 1 under the same 
conditions, except that the packet is transferred to queue 1. In a similar way 
we define the binary-valued variables D^ft), D^Ct) for the case that 
L 1 (t)>0, L 2 (t)=0 and D^(t), D^(t) for the situation when both queues are 

nonempty and departure occurs from queue 2. Then the recursive equations for 

L.(t) are : 

1 

Ljft+1) = • 


L : (t+1) = ■ 

From (Al), (A2) and (3) we have for t-*» that 

G(x,y) = F( X ,y){G(0,0) + [G(x,0)-G(0,0)][x- 1 pJ 0+ x* 1 yp} 0+ (p° Q +p} 0 )] + 

+ [G(0,y) -G(0,0) ] [y”^Q^y" lx Poi + (p 01 +P 01 )J + 

+ [G(x,y)-G(x,0)-G(0,y)+G(0,0)] [y^p^+y'^p^ + (p^+p^)]} (A3) 


Aj(t) if 

A 1 (t) + L 1 (t)-D° 0 (t)-D} 0 (t) if 
A 1 (t) + Dj 1 (t) if 

A 1 (t) + L 1 (t)+Dj 1 (t) if 

A 2 (t) if 

A 2 (t) + D} 0 (t) if 

A 2 (t)+L 2 (t)-D° 1 (t)-Dj 1 (t) if 

A 2 (t) + L 2 (t)-D° 1 (t)-Dj 1 (t) if 


LjCO-O, L 2 (t)-0 
L (t)>0, L~(t)=0 

(Al) 

L^(t)=0, L 2 (t)>0 
L 1 (t)>0, L 2 (t)>0 

Lj(t)=0, L 2 (t)«0 
(t)>0, L 2 (t)=0 

Lj^(t) =0, L 2 (t)>0 (A2) 

I-, (t)>0, L-(t)>0 


(Remember that v denotes 1-v). 
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Arranging (A5) we obtain 


Gfx v) = F(x y) b ( x >y) G ( x » Q )+ c ( x »y)G(0,y)+d(x,y)G(0,0) 
i. >yj *- x e ( x ,y) 


(A4) 


where 


b ( x >y) 

c(x,y) 

d(x,y) 

e(x,y) 


= y(p io +y p!o } - x( P?i +x Pn )+x y ( P?i + Pn-P? 0 -Pio ) 

(AS) 

= x tPorPn +x (pJrPn)] + x y(p?i + Pn-PorpJi) 

(A6) 

, o 1 , r 0 0 .1 1,. ,0 1 0 

=-y(Pio+yPio)^ x [Pii-Poi*x(Pii-p 0 i)] + x y(p 10 + p 10 + p 01 

10 1, 
+P 0l' P ll'Pll^ 

= y-F(x,y)[pJ 1 +xpJ 1 +y(pJ 1 +pj 1 )] 

(A8) 


From (A4) we see that the steady-state generating functions for the queues' lengths 
are : 

For queue 1 : 


G(x,1) 


Ffx n b(x,1)G(x,0)+c(x,1)G(0,1)+d(x,1)G(0,0) 
1 J x e(x,l) 


For queue 2 : 


G(l,y) 


FM vl b(l,y)G(l,Q)+c(l,y)GC0,y)+dCl,y)G(0,0) 
U,XJ e(l ,y) 


(A9) 


(A10) 


We still have to determine the functions G(x,0), G(0,y) and the constant G(0,0). 
To find G(0,y) let x-K) in (A3). Then 

G(0,y) = F(0,y){G(0,0) + G 1 (0,0)(p° 0+ ypJ 0 ) + 

+ [G(0,y)-G(0,0)][y- 1 p° 1+ (Po^P^)]} (All) 

where G^(x,y) is the derivative of G(x,y) with respect to x. 


From (All) we obtain 


G(0,y) 


F(0,y) 


1 - F(0,y)[y _1 p§ 1 + (pg 1 +pj 1 ) ] 


(A12) 


Before proceeding we need to prove the following Lemma. 
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Lemma 

For |x|<1 the equation 

e(x,y) = 0 

where e(x,y) is the expression of (A 8 ) has a unique solution within the unit 
circle |y|=l. 

Proof 

Since F(x,y) is not independent of x, then there exists a(i,j)>0 for some 
i and some j>0. Therefore for |x|<l and jyj = 1 we have 

|F(x,y)[p° 1 +p} 1 +y (P^-Pji)] I * 

00 00 00 00 OO CO 

s i F(x,y)| = 11 l a(i,j)x 1 y^ j $ l l aCi,j)IxJ 1 < l l (a(i,j) = l=|y| (A13) 

i=0 j=0 i=0 j=0 i=0 j=0 

Hence, applying Rouches' theorem [ 8 ], e(x,y) has exactly one zero within |y|=l for 
i x. < 1. 

Q.E.D. 

Now let t be the solution of 

FCO.tHt* 1 ?^ H- (pj^pjj)] = 1 (A14) 

in the unit circle. Using the Lemma for x=0, it is clear that t is unique. 

Then from the analyticity of G(0,y) for |v|< 1 it follows that 

,-10 0 1 
c Poi'Poi'Poi 

G (0,0)- r 1 ■■■■■ G(0,0) . (A15) 

PlO +t PlO 

Substituting (A15) in (A12), G(0,y) is determined up to the constant G(0,0). 

To find G(x,0) let f=f(x), be the solution for |x| < 1 of 

e(x,f) = 0 (A16) 

in the circle |f| =1. In the lemma it was proved that such a solution f is 


unique. 
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Then again from the analyticity of G(x,y) for |xj < 1 it follows that 


Gfx 01 - c(x,f)G(Q ,f~) +d(x,f)G(0,0) 

b(x7f) 

Substituting (A12) in (A17), the function G(x,0) 
G[0,0). To find G(0,0) let x+1 in (A9) and 
tion condition G(l,l) = 1 to obtain" : 


(A17) 

is determined up to the constant 
y-*-l in (A10) and use the normaliza- 


e 1 (l,l)-b 1 (l.l)G(l,0)+c 1 U,l)G(0,l)+d 1 (l,l)G(0,0) (A 18 ) 

and 

e 2 ^1.1) = b 2 (1,1)G(1,0)+c 2 (1,1)G(0,1)+d 2 (1,1)G(0,0) (A19) 


and from (A12) we have 


G(0,1) = F(0,1) 


p^GCO.OMp^+p^G^O.Q) 

1 - FCO.l)^ 


(A20) 


Solving (A15), (A18),(A19) and (A20), we obtain G(0,0). The ergodicity condition 
is that G(0,0) > 0. 
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Two-node networks, (a) System 1: non-symmetric ALOHA network. 

(b) System 2: tandem network; no interference at the station. 

(c) System 3: tandem network; interference at the station. 
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Fig. 4 - System 1: T versus p. 
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Fig. 14 - Comparison of approximation and simulation results of 
four-node symmetric ALOHA network. 















